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We have measured the number of Uke-sign (LS) and opposite-sign (OS) 
lepton pairs arising from double semileptonic decays of b and 6-hadrons, pair- 
produced at the Fermilab Tevatron cohider. The data samples were collected 
with the Collider Detector at Fermilab (CDF) during the 1992 - 1995 collider 
run by triggering on the existence of jifj, or efj, candidates in an event. The 
observed ratio of LS to OS dileptons leads to a measurement of the average 
time-integrated mixing probability of all produced 6-fiavored hadrons which 
decay weakly, x = 0.152 ± 0.007 (stat.) ± 0.011 (syst.), that is significantly 
larger than the world average % = 0.118 it 0.005. 
PACS number(s): 13.85.Qk, 13.20. Jf 

I. INTRODUCTION 

The time evolution of — mixing has been accurately measured in a number of 
experiments, while B^ — B^ mixing has not yet been observed. Time-independent mea- 
surements of B^ mixing offer an experimentally distinct technique to extract 5° mixing 
parameters. The time- integrated mixing probability is defined as x— r{B-^e^x) — ■> where 
the numerator includes B^ and B^ mesons and the denominator includes all B hadrons. The 
average probability is then X = fd ■ Xd + fs ■ Xs, where Xd and fd, and Xs and fs are the 
time-integrated mixing probability and the fraction of produced B^ and B^ mesons, respec- 
tively, that decay semileptonically. A measurement of x can be used to extract S° mixing 
information through Xd and Xs, or, alternatively, to extract information on the fractions of 
produced and B^ mesons. 

A precise measurement of the time-integrated mixing probability x at the Tevatron can 
also provide indications for new physics through its comparison with the LEP measurements 
and the time-dependent results from the Tevatron. For example, a recent publication [1] 
explores an explanation within the context of the minimal supersymmetric standard model 
for the long-standing discrepancy between the measured cross section for bottom-quark pro- 
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duction at the Tevatron and the next-to-leading order (NLO) prediction. Ref. [1] postulates 
the existence of a relatively hght gluino g (mass ~ 12 to 16 GeV/c^) that decays into a b 
quark and a light b squark (mass ~ 2 to 5.5 GeV/c^). The pair production of such light 
gluinos provides a bottom-quark cross section comparable in magnitude to the conventional- 
QCD component. Since ^ is a Majorana particle, its decay yields both quark and antiquark; 
therefore, gluino pair production and subsequent decay to 6-quarks will generate bb and bb 
pairs, as well as the bb final states that appear in conventional QCD production. The pair 
production of gluinos leads therefore to an increase of like-sign dileptons from weak decays 
of b quarks ^. This increase could be confused with an enhanced rate of — mixing and 
result in a value of x larger than the world average 0.118 ± 0.005 [3]. Using a previous CDF 
result [4] (x = 0.131 ± 0.020 (stat.) ± 0.016 (syst.), Ref. [1] estimates that the value of x at 
the Tevatron could be as large 0.17 ^. The x measurement in Ref. [4] is based upon muon 
pairs corresponding to an integrated luminosity of 17.4 pb~^. The present measurement, 
which makes use of a dimuon data set corresponding to an integrated luminosity of 105 pb~^ 
and an efi data set corresponding to approximately 85 pb~^, supersedes our previous result. 

In this study, the time-integrated mixing probability x is derived from the ratio of the 
observed numbers of LS and OS lepton pairs arising from 66 production. At the Tevatron, 
dilepton events result from decays of heavy quark pairs {bb and cc), the Drell-Yan process, 
charmonium and bottomonium decays, and decays of tt and K mesons. Background to 
dilepton events also comes from the misidentification of tt or X mesons. As in Ref. [4], 



^Constraints to this scenario have been derived from other data analyses (see, for example, Ref. [2] 
and experimental references therein). 

^Determinations of Xd [5], based on the direct measurement of the oscillation frequency Am^, 
are not sensitive to this type of unconventional hh production; in fact, an extra source of like- 
sign b quarks, would reduce the amplitude of the mixing asymmetry, but would not affect the 
determination of ^vrid- 
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we make use of the precision tracking provided by the CDF sihcon microvertex detector to 
evaluate the fractions of leptons due to long-hved b- and c-hadron decays, and to the other 
background contributions. 

Sections II and III describe the detector systems relevant to this analysis and the data 
selection, respectively. The analysis method, similar to the one used in Ref. [4] , is discussed 
in Sec. IV. In Sec. V, we determine the contributions of the bb and cc production to OS 
and LS dileptons. The B° — B° mixing result is derived in Sec. VI. Section VII presents 
cross-checks and studies of systematics effects. Our conclusions are summarized in Sec. VIII. 

II. CDF DETECTOR AND TRIGGER 

The CDF detector is described in detail in Ref. [6]. We review the detector components 
most relevant to this analysis. Inside the 1.4 T solenoid the silicon microvertex detector 
(SVX) [7], a vertex drift chamber (VTX), and the central tracking chamber (CTC) provide 
the tracking and momentum information for charged particles. The CTC is a cylindrical 
drift chamber containing 84 measurement layers. It covers the pseudorapidity interval [r^l < 
1.1, where r] = — ln[tan(^/2)]. In CDF, 9 is the polar angle measured from the proton 
direction, is the azimuthal angle, and r is the radius from the beam axis (z-axis). The 
SVX consists of four layers of sihcon micro-strip detectors located at radii between 2.9 and 
7.9 cm from the beam line and provides spatial measurements in the r — 4> plane with a 
resolution of 13 /xm. It gives a track impact parameter ^ resolution of about (13 -|- 40/pt) 
Hm, where pr is the track momentum measured in the plane transverse to the beam axis 
and in GeV/c units. The SVX extends ±25 cm along the z-a.xis. Since the vertex z- 
distribution for pp collision is approximately a Gaussian function with an rms width of 30 
cm, the average geometric acceptance of the SVX is about 60%. The transverse profile 

^The impact parameter is the distance of closest approach of a track to the primary event vertex 
in the transverse plane. 
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of the Tevatron beam is circular and has an rms spread of ~ 30 //m in the horizontal 
and vertical directions. The pt resolution of the combined CTC and SVX detectors is 
5ptIpt = [(0.0066)2 ^ (0.0009 (GeV/c)-^ •pr)T^^- Electromagnetic (GEM) and hadronic 
(CHA) calorimeters with projective tower geometry are located outside the solenoid and 
cover the pseudorapidity region \r]\ < 1.1, with a segmentation of Acj) — 15^^^ and Ar] — 0.11. 
A layer of proportional chambers (CES) is embedded near shower maximum in the GEM 
and provides a more precise measurement of the electromagnetic shower position. Two 
muon subsystems in the central rapidity region (|7/| < 0.6) are used for muon identification: 
the central muon chambers (GMU), located behind the GHA calorimeter, and the central 
upgrade muon chambers (GMP), located behind an additional 60 cm of steel. 

GDF uses a three-level trigger system. At the first two levels, decisions are made with 
dedicated hardware. The information available at this stage includes energy deposit in the 
GEM and GHA calorimeters, high-pr tracks found in the GTG by a fast track processor, 
and track segments found in the muon subsystems. At the third level of the trigger, events 
are selected based on a version of the off-line reconstruction programs optimized for speed. 
The lepton selection criteria used by the 3rd level trigger are similar to those described in 
the next section. 

A large fraction of the events used for this analysis are collected using two triggers that 
require two lepton candidates in an event. The first trigger requires two muon candidates; 
each muon candidate requires a track in the GTG, matched with track segments in the GMU 
system, corresponding to a particle with p^ > 2.2 GeV/c. At least one of the candidates is 
required to have track segments in both the GMU and GMP chambers. The second trigger 
requires an electron and a muon candidate. The Et threshold for the electron is 5 GeV, 
where = EsinO, and E is the energy measured in the GEM. In addition, the trigger 
requires the presence of a GTG track with px > 4.7 GeV/c and the same cf) angle of the 
GEM energy deposit. The muon candidate requires a GTG track with matched segments in 
the GMU chambers and pt > 2.7 GeV/c. 
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III. DATA SELECTION 



For this analysis we select events which contain two and only two good leptons. Good 
muons are selected by requiring Pt > 3 GeV/c and a match between the CTC track extrap- 
olated in the muon chambers and the muon segment within 3 cr in the r — (j) plane (CMU and 
CMP) and -\/l2 a in the r — z plane (CMU), where cr is a standard deviation including the 
effect of multiple scattering. In order to minimize misidentification of muons due to hadronic 
punchthrough, we require a muon segment in the CMP chambers as well as an energy de- 
posit in the calorimeters larger than 0.1 GeV but smaller than 2 and 6 GeV in the CEM 
and CHA, respectively. The identification of good electrons makes use of the information 
from calorimeters and tracking chambers. We select electrons with Et > 5 GeV, and, as in 
previous analyses [8], we require the following: (1) the ratio of hadronic to electromagnetic 
energy of the cluster, Ehad/Eem < 0.05; (2) the ratio of cluster energy to track momentum, 
E/P < 1.5; (3) a comparison of the lateral shower profile in the calorimeter cluster with that 
of test-beam electrons, Lghr < 0.2; (4) the distance between the extrapolated track-position 
and the CES measurement in the r — and z views. Ax < 1.5 cm and Az < 3.0 cm; (5) 
a comparison of the CES shower profile with those of test-beam electrons, xltrip — 
Fiducial cuts on the electromagnetic shower position as measured in the CES, are applied to 
ensure that the electron candidate is away from the calorimeter boundaries and the energy 
is well measured. Electrons from photon conversions are removed using an algorithm based 
on track information [8] . 

To ensure accurate impact parameter measurement, each lepton track is required to be 
reconstructed in the SVX with hits non-shared with other tracks in at least two layers out of 
the possible four. We also require the impact parameter of each lepton track to be less than 
0.2 cm with respect to the primary vertex ^. Lepton tracks are required to be within 5 cm 

^This cut removes most of the cosmic rays, since this background is distributed as a linear function 
of the impact parameter. 
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from the primary vertex in the z-direction. To reconstruct the primary event vertex, we first 
identify its z-position using the tracks reconstructed in the VTX detector. When projected 
back to the beam axis, these tracks determine the longitudinal position with a precision of 
about 0.2 cm. The transverse position of the primary vertex is determined for each event 
by a weighted fit of all SVX tracks which have a z coordinate within 5 cm of the vertex 
position of the primary vertex. First, all tracks are constrained to originate from a common 
vertex. The position of this vertex is constrained by the transverse beam envelope described 
above. Tracks that have impact parameter significance \d\lad^ where is the estimate of 
the uncertainty on the impact parameter d, larger than three with respect to this vertex 
are removed and the fit is repeated. This procedure is iterated until all used tracks satisfy 
the impact parameter requirement. At least five tracks must be used in the determination 
of the transverse position of the primary vertex or we use the nominal beam-line position. 
We use this procedure to avoid having the primary vertex position biased by the presence 
of heavy flavor decays [8] . The primary vertex coordinates transverse to the beam direction 
have uncertainties in the range of 10—25 /xm, depending on the number of tracks and the 
event topology. 

In the analysis, lepton pairs arising from h cascade decays are removed by selecting 
dilepton candidates with invariant mass greater than 5 GeV/c^. 

IV. METHOD OF ANALYSIS 

For leptons originating from the decay of long lived particles the impact parameter is 
d = |/37ct sin(5)|, where t is the proper decay time of the parent particle from which the 
lepton track originates, 5 is the decay angle of the lepton track with respect to the direction 
of the parent particle, and is a Lorentz boost factor. The impact parameter of the 
lepton is proportional to the lifetime of the parent particle. The markedly different impact 
parameter distributions for leptons from h decays, c decays, and other sources allow the 
determination of the parent fractions. 
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The method used to determine the bb and cc content of the data has been pioneered in 
Ref . [4] . The procedure is to fit the observed impact parameter distribution of the lepton 
pairs with the expected impact parameter distributions of leptons from various sources. 
After data selection, the main sources of reconstructed leptons are semileptonic decays of 
bottom and charmed hadrons, and prompt decays of onia and Drell-Yan production. 

Monte Carlo simulations are used to model the impact parameter distributions for leptons 
from b and c decays. We use the HERWIG Monte Carlo generator program [9] to generate 
hadrons with heavy flavors ^, the QQ Monte Carlo program [10] to decay hadrons with heavy 
flavor, and the QFL Monte Carlo simulation of CDF [8] to model the detector's response. 
Impact parameter distributions for simulated b and c decays are shown in Figures 1(a) 
and (b), respectively. Since lifetimes of bottom and charmed hadrons (ctb — 480 /xm and 
ctd — 200 /xm) are much larger than the average SVX impact parameter resolution in these 
data sets (~ 15 //m), the dominant factor determining the impact parameter distribution 
is the kinematics of the semileptonic decays which is well modeled by the simulation (see 
Sect. VII). The fraction of leptons from sequential b decays (6 — > cX,c — > lY) is also 
determined with the simulation. Leptons from sequential b decays have slightly different 
kinematics and slightly larger ct than leptons coming from direct b decays; these two effects 
compensate and the simulated impact parameter distribution of leptons from sequential 
decays is indistinguishable from that of leptons from direct b decays. 

The impact parameter distribution of leptons from prompt sources such as quarkonia 
decays and Drell-Yan production is plotted in Fig. 1(c) and is derived using muons from 

^We use option 1500 of version 5.6, generic 2—^2 hard scattering with > 5 GeV/c, with the 
same setting of the herwig parameters used in Ref. [8] . In the generic hard parton scattering, bb 
and cc pairs are generated by HERWIG through processes of order (LO) such as gg bb (direct 
production). Processes of order are implemented in HERWIG through flavor excitation processes, 
such as gb — > gb, or gluon sphtting, in which the process gg — > gg is followed hy g ^ bb. 
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T(1S) decays ^ (see Fig. 2). 

Lepton tracks from tt and K in-flight decays are also regarded as prompt tracks since 
the track reconstruction algorithm rejects tracks with appreciable kinks. Tracks of tt and K 
mesons, which mimic the lepton signal, are also regarded as prompt since the average heavy 
flavor contribution per event is neghgible (see Sect. VII). 



^We use templates derived from the data to account properly for non-Gaussian tails of the impact 
parameter distribution. The impact parameter distribution of electrons from a smaller statistics 
sample of Z e+e~ is also well modeled by the muon template. 
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FIG. 1. Impact parameter distributions of leptons coming from b decays (a), c decays (b), and 
prompt leptons (c). Distributions are normalized to unit area; differences between n — ji and 
e — jJL templates are due to the different pT thresholds. The ratio of the number of events with 
d < 0.008 cm to that with d > 0.008 cm is 1.04, 2.85, and 32.3 for the histograms (a), (b), and (c), 
respectively. 
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FIG. 2. Invariant mass distribution of OS dimuons in the T region. The impact parameter 
distribution in Figure 1(c) is derived using muons with invariant mass between 9.28 and 9.6 GeV/c^. 
The background is removed using dimuons with invariant mass between 9.04 and 9.2 GeV/c^ and 
between 9.54 and 9.7 GeV/c^. Dimuon events in the mass range 9.2 — 10.5 GeV/c^, which are 
dominated by T production, are not used in the x analysis. 



Since there are two leptons in an event, the fit is performed in the two-dimensional space 
of impact parameters. Each axis represents the impact parameter of one of the two leptons. 
In filling the histograms, the lepton ordering by flavor type or transverse momentum is 
randomized. The two-dimensional impact parameter technique exploits the fact that the 
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lepton impact parameters are independent uncorrelated variables ^. The two-dimensional 
template distributions for each type of event are made by combining the relevant one- 
dimensional distributions in Fig. 1. 

A binned maximum log likelihood method is used to fit simultaneously the impact pa- 
rameter distributions of OS and LS dileptons. The likehhood L is defined as 



where n{i,j) is the number of events in the {i,j)th bin. The function l^j is defined as 



kj = BB ■ S,{i) ■ S,{j) + CC ■ Sc{i) ■ Sc{j) + PP ■ S^{i) ■ Sp{j) + 

0.5 • [BP ■ {S,{{) ■ S,{j) + S,{t) ■ S,{j)) + CP ■ • S,{j) + S,{t) ■ S,{j))] 

where Sh, Sc, and Sp are the impact parameter templates shown in Fig. 1(a), (b), and (c), 
respectively. The fit parameters BB, CC, and PP represent the bb, cc and prompt dilepton 
contributions, respectively. The fit parameter BP (CP) estimates the number of events 
in which there is only one b (c) quark in the detector acceptance and the second lepton is 
produced by the decay or the misidentification of tt and K mesons ^. Figure 3 compares 
projections of the two-dimensional distributions for each type of dilepton contribution to 
the likelihood. Because of sequential decay and mixing, the bb production results in both 
OS and LS dileptons. For LS dileptons, one expects no contribution from cc production. 

We do not fit dimuon events with invariant mass between 9.2 and 10.5 GeV/c^ since OS 
dimuons are dominated by T meson production. The PP contribution to e/i events can 



, IS approx- 



^The correlation between the two impact parameters, p = //("^i <di>){d2^ <d2>)SdiSd2 ^ 
imately 0.04 in the data samples and their heavy flavor simulations. 

^According to the simulation, supported by the measurement in Ref. [11], approximately 90% of 
the bb and cc events with an identified lepton from heavy flavor decay do not contain the second 
heavy flavored hadron in the detector acceptance. Therefore, we ignore the small contribution to 
misidentified leptons due to tt and K mesons from heavy flavor decays (see Sec. VII). 
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only arise from misidentified leptons (rr Drell-Yan production is negligible) and is expected 
to be equal for OS and SS dileptons. Therefore, in the fit to e/j, data, the PP components 
in OS and LS dileptons are constrained to be equal within the statistical error (technically, 
we add the term 0.5 x (PP(OS) - PP{LS)f /{{PP{OS) + PP(LS)) to the function - InL 
used by the fit). In dimuon events, where the Drell-Yan contribution is relevant, OS leptons 
have a larger PP component than LS dileptons. The BP and CP contributions, in which 
one lepton is fake, are expected to be the same for OS and LS dileptons, and in the fit 
are constrained to be equal within the statistical error. One also expects the BP and CP 
contributions to have approximately the same size ^. 



^According to the simulation, the cross section for producing at least one c hadron in the detector 
acceptance is approximately a factor of two larger than the cross section for producing at least 
one b hadron in the detector acceptance. Since the efficiency for detecting a lepton from a c decay 
is approximately 40% of that for detecting a lepton from a b decay, one expects the bb and cc 
contributions to events with at least one identified lepton to be approximately equal. In contrast, 
the bb and cc cross sections for producing events which contain 2 hadrons with heavy flavor in the 
detector acceptance are dominated by the LO term and are approximately equal; one therefore 
expects the bb contribution to dilepton events to be much larger than the cc contribution. 



18 




0.1 0.15 
d (cm) 



0.1 0.15 
d (cm) 



0.1 0.15 
d (cm) 



0.1 0.15 
d (cm) 



0.2 




0.2 



FIG. 3. Projections of the two-dimensional impact parameter distributions of the different com- 
ponents used to fit the dimuon data (see text). The top- left distribution shows the shapes of the 
prompt, b and c templates used to construct the different two-dimensional distributions used in 
the likelihood function. All distributions are normalized to unit area. 
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V. RESULT 



We show the result of the fit to the data for dimuon and eji events in subsections A and 
B, respectively. 

A. Dimuon events 

The observed two-dimensional impact parameter distributions for OS and LS dimuons 
arc plotted in Figure 4. We do not use dimuon events with invariant mass between 9.2 and 
10.5 GeV/c^ since OS are largely dominated by T meson production. There are 18420 OS 
dimuons and 9279 LS dimuons after the removal of 6264 OS and 1302 LS dimuons with 
invariant mass in the T region. 
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FIG. 4. Two dimensional impact parameter distributions for (a) OS and (b) LS dimuons. 



One sees that a handful of events in Fig. 4(a) cluster along the diagonal line di — 6,2. 
These events are due to cosmic rays. We minimize their contribution by fitting only events 
with di + d2 < 0.2 cm. As shown in Sec. VII, the fit result is unaffected by the inclusion 

20 



of events with di + d2 > 0.2 cm. When all the hkelihood terms are used to fit the data, 
the best fit, as expected, returns CC = ± 40 LS events. However, while the fit finds 
an appreciable BP component, it returns CP = ± 110 in both LS and OS events. When 
fitting the data with all components, the fit gets blocked when hmiting the CC{LS) and CP 
parameters to positive values, and it returns reliable errors only when allowing the CC and 
CP terms to have also unphysical (negative) values. Since these unphysical values produce 
an overestimate of the size and the error of the remaining components, we fit again the data 
setting to zero the CC term in LS events and the CP contribution to OS and LS events 

The fit result is shown in Table I. The parameter correlation matrix is listed in Table II. 
The best fit returns — InL = 3076. The probability of the — InL value returned by the fit 
is determined by fitting Monte Carlo pseudoexperiments. In each experiment, we randomly 
generate difi^erent components with average size as determined by the fit to the data and 
allowing for Poisson fiuctuations; the impact parameter distribution for each component is 
randomly generated from the corresponding templates used in the fit. We find that 40% 
of the fits to the pseudoexperiments return a — In L value equal or larger than 3076. For 
a comparison of the data and the fit results, projections of the two-dimensional impact 
parameter distributions are shown in Fig. 5. Since the fit appears to underestimate the data 
for di > 0.12 cm, we have fitted the data excluding points at impact parameters larger than 
0.12 cm; this fit returns a result identical to that of the standard fit. 

Using Table I, one derives a ratio of LS to OS dimuons due to bb production which is 
i? = 0.537 ±0.018. 



^°In Sec VII, we show that this happens in 15% of simulated pseudo-experiments due to the fact 
that CC, BP and CP templates are quite similar. In addition, we show that the fit result does 
not vary when constraining the BP and CP components to be, as expected, equal within their 
statistical error. 
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FIG. 5. The projection of the impact parameter distribution of (a) OS and (b) LS dimuons onto 
one of the two axis is compared to the fit. 

TABLE I. Number of events attributed to the different sources of dimuons by the fit to OS and 
LS dimuons with di + d2 < 0.2 cm. The errors correspond to a 0.5 change of the — InL. 



Component OS LS 

BB 10476 ± 223 5630 ± 132 

CC 2469 ± 360 

PP 3603 ± 161 1914 ± 87 

BP 1566 ± 165 1555 ± 157 

CP 
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TABLE II. Parameter correlation coefHcients returned by the fit listed in Table I. 

Component BB{OS) CC{OS) PP{OS) BP{OS) BB{LS) PP{LS) 

CC{OS) -0.70 

PP{OS) 0.53 -0.73 

BP{OS) -0.03 -0.46 0.05 

BB{LS) 0.02 0.31 -0.03 -0.66 

PP{LS) 0.02 0.27 -0.03 -0.58 0.25 

BP{LS) -0.03 -0.44 0.05 0.94 -0.71 -0.62 
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B. e/x events 



Figure 6 shows the observed two-dimensional impact parameter distributions for OS and 
LS eii pairs. There are 7802 OS and 4331 LS e/x events 
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FIG. 6. Two dimensional impact parameter distributions for (a) OS and (b) LS eji events. 



When all the likehhood terms are used to fit the data, the best fit, as expected, returns 
CC = ± 80 LS events. However, while the fit finds an appreciable BP component, it 
returns CP = ± 130 in both LS and OS events. As in the case of dimuon events, the 
fit gets blocked at the lower limits when the CC{LS) and CP parameters are bound to 
be positive, and we exclude these terms in the fit hkelihood. The fit result is shown in 
Table III and the parameter correlation matrix is listed in Table IV. The best fit returns 



Since lepton tracks are reconstructed requiring at least two hits in the SVX detector close to the 
beam pipe, the number of electrons due to unidentified photon conversion is negligible (no larger 
than three). 
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— InL = 2481. As for dimuon events, the probability of the — InL value returned by the fit 
is determined by fitting Monte Carlo pseudoexperiments. We find that 62% of the fits to the 
pseudoexperiments return — InL values equal or larger than 2481. For a comparison of the 
data and the fit result, projections of the two-dimensional impact parameter distributions 
are shown in Fig. 7. Since the fit appears to underestimate the data for di > 0.1 cm, we have 
fitted the data excluding points at impact parameters larger than 0.1 cm; this fit returns a 
result identical to that of the standard fit. 

Using Table III one derives that the ratio of LS to OS dileptons due to bb production is 
R = 0.560 ± 0.024. 
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FIG. 7. The projection of the impact parameter distribution of (a) OS and (b) LS e// pairs onto 
one of the two axis is compared to the fit. 
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TABLE III. Number of events attributed to the different sources by the fit to OS and LS e/i 

pairs. The errors correspond lo a 0.5 cliauge of I lie — luZ. 



Component OS LS 

BB 5099 ± 138 2852 ± 90 

CC 1126 ± 162 

PP 906 ± 60 875 ± 52 

BP 536 ± 107 529 ± 102 

CP 



TABLE IV. Parameter correlation coefficient returned by the fit listed in Table III. 

Component BB{OS) CC{OS) PP{OS) BP{OS) BB{LS) PP{LS) 

CC{OS) -0.63 

PP{OS) 0.38 -0.37 

BP{OS) -0.23 -0.33 -0.43 

BB{LS) 0.12 0.29 0.18 -0.67 

PP{LS) 0.31 -0.14 0.76 -0.56 0.23 

BP{LS) -0.23 -0.29 -0.45 0.95 -0.70 -0.59 



26 



VI. AVERAGE MIXING PROBABILITY 



The average B^B^ mixing probability is defined as 

r(^o ^B^^ l+X) 

^ ~ r{B i^x) 

where the numerator includes B^ and B^ mesons and the denominator includes all B 
hadrons. In absence of mixing, the double semileptonic decay of a BB pair results in 
an OS lepton pair; when one of the mesons undergoes mixing a LS lepton pair is produced. 
The mixing probability x can therefore be inferred from i?, the ratio of LS to OS dileptons 
due to bb production. 

The sequential decays of 6-hadrons also contribute to R. The fraction of leptons from 
sequential decays, is evaluated using the simulation. Using simulated dimuon events, we 
find = 0.123 with a 12% uncertainty As for the study of Ref. [4], the uncertainty on 

comes from the uncertainty of the relative branching ratios of h and c semileptonic decays 
(±11%) and the uncertainty of the detector acceptance for sequential leptons with respect 
to that for leptons from direct decays (±6%). Using the e/x simulation, we derive /e = 0.060 
and = 0.142 with a ±12% systematic uncertainty. 

The ratio R is related to the time-integrated mixing probability in the following way: 

P_ /[X^±(1-X)^]±2x(l-X)(l-/) 
(l-/)[x2±(l-x)^]±2x(l-x)/ 

where / = 2/^(1 -/^) = 0.2157± 0.0226 (syst.) for dimuon events and / = /e±/^-2/e/^ = 
0.1850 ± 0.0204 (syst.) for eji events. Systematic errors due to other sources are neghgible 
with respect to that arising from the / uncertainty, and are neglected (see Sec. VII). 
Prom the observed values of i?, we derive the following mixing probabilities: 
X = 0.136 ± 0.009 (stat.) ± 0.014 (syst.) for dimuon events 
X = 0.165 ± 0.011 (stat.) ± 0.011 (syst.) for e/x events 



^^Technically this fraction accounts also for the 0.4% fraction of events which contain more than 
two hadrons with heavy flavor. 
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Since we use events containing two and only two leptons, the results from the dimuon and 
e/x data sets are statistically independent. Therefore, we combine the two results and derive 
an average mixing probability x = 0.152 ± 0.007 (stat.) ± 0.011 (syst.) 

This value of the mixing probability agrees with all previous results from pp colliders 

X = 0.157 ± 0.020 (stat.) ± 0.032 (syst.) (UAl [12]) 

X = 0.176 ± 0.031 (stat. + syst.) ± 0.032 (model) (CDF [13]) 

X = 0.131 ± 0.020 (stat.) ± 0.016 (syst.) (CDF [4]) 
but is significantly larger than the world average x = 0.118 ± 0.005 [3], which is dominated 
by the LEP measurements at the Z-pole Since our result is statistically very different 
from the world average, we have investigated the error behaviour beyond one a. For an 8 
unit increase of the — InL value (4 a uncertainty), the errors of the BB{OS) and BB{LS) 
terms returned by the fit increase by a factor of four, and we derive a 4 cr statistical error 
of 0.029 for the combined value of x- 

VII. CROSS-CHECKS OF THE RESULT AND STUDY OF ADDITIONAL 

SYSTEMATIC EFFECTS 

In this section, we first perform several cross-checks of the x result, and then investi- 
gate its sensitivity to the modeling of the production and weak decay of heavy quarks. In 
subsection A we verify that the ratio of the number of lepton pairs due to cc production 
to that due to bb production returned by the various fits is consistent with the theoretical 
expectation. Subsection B compares our result to the previous CDF measurement, which 
used a subset of the data available for this analysis. Subsection B also verifies that the x 



^^The systematic error quoted in Ref. [4] (±0.016) is larger to account for the fact that the BP 
and CC terms are not fitted independently. 

^^The world average assumes that the fractions fd and fg at the Tevatron are equal to those at 
the Z-pole. 
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result is not affected by the small cosmic ray background present in the dimuon data sample. 
Subsection C shows that the x result is not affected by the fact that we have excluded the 
CP component in the fit likelihood. Subsections D, E and F explore the dependence of our 
result on the mixture of the different b and c hadrons, on the ratio of bb to cc production 
cross section, and on the transverse momentum distribution of hadrons with heavy flavor 
predicted by the QCD simulation. In analogous measurements, these effects are usually 
not considered since they are hard to quantify and to implement consistently into the QCD 
generator. We investigate them either by changing the heavy flavor composition of the data 
with proper kinematical selections, or with reasonable modifications of the simulation pre- 
diction. Finally, subsections G and H verify the templates used to separate the contribution 
of semileptonic decays of heavy flavor from that of leptons due to misidentified hadrons or 
prompt sources as the Drell-Yan process. We show that all above effects change our result 
by a very small fraction of the quoted statistical and systematic errors. We report changes 
in R when the sequential fraction fi is not affected by the particular study, and also changes 
in X when /; is affected; a summary of the different results is presented in subsection I. 



A. Ratio of the cc to bb production 

The difference between the x measurements at the Tevatron and LEP may not require 
an explanation in terms of new physics; however, if we entertain the hypothesis [1] that the 
enhancement of the bb cross section at the Tevatron with respect to the NLO prediction 
may be caused by pair production of light gluinos decaying to a bottom quark and a bottom 
squark, which in turn produces an apparent increase of x with respect to LEP, then the 
ratio of the cc to bb cross sections should be approximately a factor of two smaller than 
what is predicted by the Standard Model. Therefore, it is of interest to compare the ratio 
of the numbers of leptons due to cc and bb production in the data and the simulation. 

The dimuon fit in Table I returns a ratio ^ = 0.15 ± 0.02 (stat.). In the simulation, 
this ratio is 0.18 ± 0.02 (stat.). 
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The fit to efj, data in Table III returns a ratio ^ = 0.14 ±0.02 (stat.). In the simulation, 
the ratio is 0.12 with a negligible statistical error. 

As shown in Ref. [14], which studies events with jets corresponding to partons with 
transverse momentum larger than 20 GeV/c, the HERWIG generator predicts heavy flavor 
cross sections which are approximately a factor of two larger than the NLO calculation [15] 
and models correctly the cc and bb cross section observed at the Tevatron. However, muons 
in the present analysis correspond to partons with pt > 6.5 GeV/c (electrons to partons 
with Pt > 9 GeV/c). A priori, there is no guarantee that HERWiG still does a good job 
in predicting the ratio ^ also in this data set which corresponds to a hard scattering 
with smaller transverse momenta (the inclusive bb and cc cross sections are approximately 
a factor of 40 larger in this data set than in the jet data studied in Ref. [14]). We cross- 
check the ratio of the cc to bb parton-level cross sections evaluated with HERWiG with two 
different NLO Monte Carlo calculations. In HERWIG, the ratio of the cc to bb cross sections 
for producing both heavy quarks with [r/l < 1 and transverse momentum large enough to 
produce an electron with Et > 5 GeV and a muon with pt > 3 GeV/c is 1.37. In the mnr 
calculation [15], this ratio is found to be 1.39, while the CASCADE Monte Carlo generator [16] 
predicts a value of 1.39 [17]. 

We conclude that the ratio of dileptons due to cc production to that due to bb production 
at the Tevatron is consistent with the prediction of the presently available Monte Carlo 
generators. 

B. Cosmic ray background in dimuon events and comparison with the previous CDF 

result 

The previous CDF measurement of x [4] uses a subset (17.4 pb~^) of the dimuon sample 
(105 pb~^) collected by CDF and used in the present analysis. There are minor differences 
in the data selection. In the present analysis we exclude dimuons with impact parameters 
di + d2 > 0.2 cm to reduce the impact of the cosmic ray background, and we exclude the T 



30 



invariant mass region which has a neghgible fraction of heavy flavor contribution. 

To study our sensitivity to the cosmic ray background we have performed a fit to the data 
which includes dimuons with di+d2 > 0.2 cm. This fits returns a ratio R = 0.533±0.018 (the 
standard fit yields R — 0.537 ± 0.018). We conclude that the small cosmic ray background 
does not affect the fit result. 

In order to compare with the result in Ref. [4] we fit the data including the T mass 
region. Because of the slightly different selection, the total number of events in the present 
analysis, 35265, is 24% larger than the number of events selected in Ref. [4] (4750 events) 
multiplied by the ratio of the relative luminosities. The fit which includes this mass region 
is shown in Table V. The fit returns a total of 18737 ± 275 dimuon events due to bb 
production. Consistently, this number is 25% larger than the number of dimuon events 
attributed in Ref. [4] to bb production (2471 ± 104 events) multiplied by the ratio of the 
relative luminosities. This fit that includes the T mass region yields R — 0.535±0.017 (stat.), 
which compares well to the result of our standard fit and the value R — 0.502 ±0.041 (stat.) 



TABLE V. Number of events attributed to the different sources of dimuons by the fit to OS and 
LS dimuons including the invariant mass region between 9.2 and 10.5 GeV/c^. 



in Ref. [4]. 



Component 



OS 



LS 



BB 



12202 ± 237 



6535 ± 139 



CC 



2849 ± 388 







PP 



7601 ± 189 



2173 ± 94 



BP 



1662 ± 175 



1658 ± 167 



CP 
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C. Effect of neglecting the CP component in the likelihood function 



In order to estimate correctly the uncertainties of the bb and cc contributions returned 
by the fit, we had to set to zero the CP component, which is expected to be of the same 
size of the BP component ^. We have performed a number of pseudo-experiments of ap- 
proximately the same size and composition as the data. In each pseudo-experiment, the 
impact parameters of the dileptons contributed by a given component are extracted from 
the corresponding two-dimensional template used to fit the data. Each pseudo-experiment 
has been fitted as the data, and the result of 125 pseudo-experiments is shown in Table VI. 
In 15% of the pseudo-experiments, the CP value returned by the fit is so close to zero that 
the fit gets blocked at the lower limit; as for the data, the CP term has to be ignored in the 
likelihood in order to estimate correctly the uncertainty of the BB term. 

TABLE VI. Number of generated and fitted events in 125 pseudo-experiments. We list the 
average and the rms spread of the values returned by the fits. 

Component generated fitted 



BB 



8000 



7998 ± 247 



CC 



4000 



3991 ± 544 



PP 



4000 



3999 ± 348 



BP 



1200 



1204 ± 505 



CP 



1200 



1196 ± 812 
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We have further investigated the sensitivity of the R result to the value of the CP 
component returned by the fit by constraining it to be equal to the BP contribution within 
the statistical error. The fit results are shown in Table VII for dimuon events and in 
Table VIII for efj, events. These fits return R — 0.533 ± 0.016 (the standard fit returns 
R = 0.537 ± 0.018) for dimuon events and R = 0.559 ± 0.023 (the standard fit returns 
R = 0.560 ± 0.024) for e/x events. 

TABLE VII. Number of events attributed to the different sources of dimuons by the fit to OS 
and LS dimuons with di + d2 < 0.2 cm. The errors correspond to a 0.5 change of — InL. 
Component OS LS 

BB 10691 ± 232 5695 ± 134 

CC 2203 ± 404 

PP 3328 ± 166 1536 ± 122 

BP 1009 ± 130 1001 ± 126 

CP 878 lb 122 869 ± 117 



TABLE VIII. Number of events attributed to the different sources by the fit to OS and LS e/i 

events. The errors correHpoud lo a 0.5 ehrUige of —hiL. 

Component OS LS 

BB 5171 ± 134 2892 ± 92 

CC 1083 ± 162 

PP 798 ± 70 767 ± 64 

BP 312 ± 63 308 ± 60 

CP 300 ± 61 293 ± 58 
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D. Sensitivity to the b and c lifetime 



The impact parameter distribution of leptons from b and c decays has some dependence 
on the hfetime uncertainty. We have varied the average 6-hadron hfetime in the simulation 
by ±10% and refit the data with the resulting templates in order to investigate which effect 
might have the possibility that the relative fractions of different 6-hadrons in the simulation 
are grossly different from the data. The fractions of the BB components, which are returned 
by the fit, change by approximately ±9% for both OS and LS dileptons; however, the ratio 
R changes by less than 0.2%. 

Since cc events contribute only to OS events, we have studied the sensitivity of the fit 
to the impact parameter template for c semileptonic decays. We have constructed impact 
parameter templates by varying in the simulation the relative ratio of D"^ to mesons by 
±30% The CC component in OS dileptons returned by the fit changes by approximately 
±10%. In the fit, this change is mostly compensated by the BP component, and the BB 
contribution to OS dilepton changes by less than ±0.1%. 

E. Sensitivity to the cc contribution 

The cc production contributes only OS dileptons. The value of R returned by the fit 
can be affected by a poor modeling of this contribution. We investigate this possibility by 
analyzing a data sample with a smaller fraction of cc contribution. According to the HERWIG 
generator program, and also to the mnr Monte Carlo program [15], the ratio of the cc to bb 
cross sections for producing both heavy flavor partons with l^y] < 1 and transverse momenta 
larger than 9 GeV/c is 1 while in the simulation of the standard e/j, data set is 1.37. 

This kinematical situation is modeled by selecting muons, as well as electrons, with 
> 5 GeV/c. We derive from the simulation of this data set new impact parameter 



The lifetime is cr = 315 iim for the meson and cr = 123 fim for the meson. 
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templates for b- and c-hadron decays. The fit result is shown in Table IX. The fit yields 

R — 0.524± 0.034. In this case, the fractions of sequential decays are /e = 0.060, = 0.092, 

and / = 0.1410 ± 0.0158 (syst.). It follows that x = 0-170 ± 0.015 (stat.) ± 0.007 (syst.), in 

agreement with the result of the standard fit x = 0.165 ± 0.011 (stat.) ± 0.011(syst.). 

TABLE IX. Number of events attributed to the different sources by the fit to OS and LS e/i 

events in which both leptons have pt > 5 GeV/c. The errors correspond to a 0.5 change of — InL. 

Component OS LS 

BB 2113 ±86 1107 ±57 

CC 421 ±98 

PP 265 ± 36 249 ± 31 

BP 163 ± 68 159 ± 65 

CP 



F. Sensitivity to the modeling of the kinematics 

Because we select leptons above a certain pt threshold, the impact parameter templates 
for leptons from semileptonic decays of heavy flavors have some dependence on the modeling 
of the Pt distribution of the parent hadron with heavy flavor The modeling of the pt 
distribution of the parent hadron with heavy flavor can be affected by a wrong estimate of 
the relative contribution of processes of order and a^, or by an incorrect modeling of the 
hadronization of heavy quarks In the next two subsections, we investigate the sensitivity 
of our result to these effects. 



^^In the extreme case of a lepton with pT close to the 5 GeV/c threshold, parent hadrons with a 
5 GeV transverse energy produce leptons with zero impact parameter. 

^^In the simulation partons arising from diagrams are slightly stiffer than those contributed 
by diagrams. 
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1. Dileptons with 6(f) > 2.4 



According to the simulation, the fractional contribution of bb and cc direct production 
(LO term) increases with increasing S(f), the azimuthal opening angle between the two lep- 
tons. Using dileptons with d4> > 2.4 rad, the number of simulated events due to bb and cc 
production is reduced by 64% and 66%, respectively. At the same time, the fraction of direct 
production in bb events increases from 71% to 84% and the fraction of direct production in 
cc events increases from 66% to 76%. 

Using this selection, the data consist of 4872 OS and 2745 LS dileptons. The result of the 
fit to these events using standard templates is shown in Table X. We derive R — 0.576±0.032, 
in good agreement with the standard fit result R = 0.560 ± 0.024. 

TABLE X. Number of events attributed to the different sources by the fit to OS and LS e/i 
events with Scf) > 2.4. The errors correspond to a 0.5 change of — In L. 

Component OS LS 



BB 



3255 ± 110 



1874 ± 75 



CC 



688 ± 129 







PP 



534 ± 47 



513 ±41 



BP 



314 ± 88 



310 ± 84 



CP 
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2. Dependence on the pr spectrum of the parent hadron with heavy flavor 

As shown by Fig. 21 of Ref. [18] and Figures 7 and 8 of Ref. [14], our simulation models 
quite well the hadronization of b and c quarks with transverse energy larger than 20 GeV. 
As shown in Fig. 8, the simulation also models correctly the lepton transverse momentum 
distributions in the e/i data. Because the lepton distribution depends on the px distribution 
of the parent parton and its fragmentation function, we use a comparison between data 
and simulation to evaluate their global uncertainty. A fit of the lepton pt spectra with the 
simulated shapes weighted with the function pj., where o; is a free fit parameter, returns 
a = 0.003 ± 0.023. In the simulation, such changes of lepton px distributions can be 
modeled by reweighting the pr distribution of the parent parton with the function p^^~^'^. 
Fits to the e/i data using templates constructed with these modified simulations return 
R = 0.557 ± 0.024 for P = 0.05 and R = 0.559 ± 0.024 for /3 = -0.05 (the result of the 
standard fit is 0.560 ± 0.023). 
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FIG. 8. Comparison of the transverse momentum distributions of electrons (a) and muons (b) in 
the data and in the heavy flavor simulation. The bottom plot (c) shows the transverse momentum 
distribution of all other tracks in e\x events. Data and simulation are normalized to the same 
number of events. 
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G. Dependence on the modeling of the impact parameter distributions 

For tracks in a jet, the impact parameter resolution in the data is shghtly larger than in 
the parametrized QFL detector simulation which has in input the SVX-hit resolution of the 
data [8] . This is believed to be due to the probability of reconstructing a track with spurious 
SVX hits, which in the data is larger than in the simulation because the SVX occupancy in 
the data is also larger. In JET 20 data the transverse energy deposited by charged tracks 
in a cone of radius 0.2 in the rj — (p space around the axis of a lepton contained in a jet is 
~ 18 GeV. For the events used in this analysis, the transverse energy deposited by charged 
tracks in a cone of radius 0.2 around each lepton is ~ 0.8 GeV; in this case, the transverse 
momentum distribution of all charged tracks in the dilepton events, plotted in Fig. 8(c), is 
also well modeled by the simulation. 

To further investigate the sensitivity to spurious SVX hits, we have repeated our study 
by using only leptons with 4 SVX hits; we also require that at least two of the hits are 
not shared with other tracks. We also make use of new templates for prompt leptons, and 
leptons from b- and c-hadron decays constructed using this track selection. 

With this selection, the dimuon data consist of 9822 OS and 4785 SS pairs. Table XI lists 
the result of the fit to dimuon events passing this selection. The fit yields R — 0.548 ±0.025, 
in good agreement with the result of the standard fit i? = 0.537 ± 0.018. 

The e/j, data consist of 4465 OS and 2355 SS pairs with 4 SVX hits. Table XII lists the fit 
result. The fit yields R — 0.559 ± 0.029, in good agreement with the result of the standard 
fit i? = 0.560 ± 0.024. For a comparison of the data and the fit results, projections of the 
two-dimensional impact parameter distributions are shown in Fig. 9. The combined result 
yields an average mixing parameter x — 0.154±0.009 (stat.) ±0.011 (syst.), to be compared 
to the standard fit result x = 0.152 ± 0.007 (stat.) ± 0.011 (syst.). 



^^Events collected with a trigger that requires at least one jet with Et >20 GeV. 
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TABLE XI. Number of events attributed to the different sources by the fit to OS and LS dimuons 

with i SVX hilti and di + (I2 < 0.2 cm. The errors correspond to a 0.5 change (^f ihe — hiZ. 

Component OS LS 

BB 4990 ± 150 2735 ± 90 

CC 1818 ± 245 

PP 2237 ± 112 1289 ±63 

BP 740 ± 110 743 ±106 

CP 



TABLE XII. Number of events attributed to the different sources by the fit to OS and LS e/x 

events \vi(h 1 SVX hits. Th(^ (ni'ors corr(^si")()nd lo a 0.5 change of lh(^ — Ini. 

Component OS LS 

BB 2768 ± 99 1547 ± 66 

CC 831 ± 121 

PP 575 ± 44 552 ± 37 

BP 266 ± 76 264 ± 73 

CP 
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FIG. 9. The projection of the impact parameter distributions in the data (•) is compared to the 
the fit result (dashed histograms). 



H. Leptons faked by tracks from hadronic decays of hadrons with heavy flavor 

In the standard fit to the data, we have approximated the impact parameter distribution 
of fake leptons with that of leptons from prompt sources. The fits return a BP component 
which is 15% (dimuon events) and 10% (e/x events) of the BB component. According to the 
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simulation, only 7.5% of the events due to the BP component contain a second hadron with 
heavy flavor which decays hadronically; in these events, less than 50% of the tracks, which 
are fake-lepton candidates, arise from the decay of the heavy flavored hadron; in addition, 
80% of the lepton faked by tracks from hadronic decays of heavy flavors carry a charge with 
the same sign of that of the parent heavy flavor quark. Therefore, one estimates that the 
effect of this approximation on R is of the order of 10~^ 

We cross-check our conclusion by modehng fake leptons with new templates, called F 
(instead of P) , derived in a sample with a comparable contamination of hadrons with heavy 
flavor. This sample consists of events containing a jet with Et > 20 GeV. As shown by the 
study in Ref. [8], JET 20 data contain a 9.5% fraction of heavy flavor. After removing events 
in which jets contain a soft lepton (SLT tag) or a displaced secondary vertex (SECVTX 
tag), the contamination of heavy flavor is 7.1% (comparable to the fraction of heavy flavor 
with hadronic decay contributing to the BF and CF components). The new template is 
constructed by using all tracks with Pt > 3 GeV/c and pointing to the CMUP flducial 
volume. Figure 10 compares the new template to the one derived using prompt muons. 



^This is supported by the fact the CC component in LS dilepton, which can only be contributed 
by leptons faked by tracks from hadronic decays of charmed hadrons, is found negligible by our fit 
with a 1 cr upper limit of 1.6% of the CC contribution to OS dileptons. 
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FIG. 10. Comparison of the impact parameter distributions of lepton candidate tracks in JET 20 
data and of leptons coming from T(15) decays. 

Tables XIII and XIV list the results of the fits to dilepton events with 4 SVX hits 
when using templates which account for the heavy flavor contribution to fake leptons. The 
fits return R = 0.570 ± 0.027 for dimuon events, and R = 0.562 ± 0.034 for e/i events. The 
combined result yields an average mixing probability x = 0.159±0.010 (stat.) ±0.011 (syst.) 
to be compared to the standard fit result x = 0.154 ± 0.009 (stat.) ± 0.011 (syst.). 
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TABLE XIII. Number of events attributed to the different sources of dimuons by the fit to 
OS and LS dimuons with 4 SVX hits and di + ^2 < 0.2 cm. Fake leptons for the BF and CF 

components are modeled with a template derived in JET 20 data. 

Component OS LS 

BB 4781 ± 150 2723 ± 90 

CC 2207 ± 222 

PP 2018 ± 111 1251 ± 64 

BF 787 ± 108 796 ± 104 

CF 



TABLE XIV. Number of events attributed to the different sources of dimuons by the fit to OS 
and LS e/x with 4 SVX hits. Fake leptons for the BF and CF components are modeled with a 
template derived in JET 20 data. 

Component OS LS 

BB 2743 ± 103 1541 ± 68 

CC 857 ±118 

PP 586 ± 45 566 ± 38 

BF 257 ± 76 256 ± 73 

CF 
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I. Summary of the cross-checks 



Table XV lists the x values resulting from the different cross-checks presented in this 
section. All x nieasurements are consistent with the main result presented in Sec. VI. 



TABLE XV. Summary of the cross-checks presented in Sec. VII. 


The X error is statistical only. 


Data set 


fit type 


X 


H/j, + e/j, 


standard 


0.152 ± 0.007 


lifj, + e/j, 


BP = CP (Sec. VII C) 


0.151 ±0.007 


fi/j, + e/j, 


4 SVX hits (Sec. VUG) 


0.154 ±0.009 


li/j, + e/j. 


4 SVX hits, JET 20 fakes (Sec. VII H) 


0.159 ±0.010 


eji 


standard 


0.165 ±0.011 


efi 


A(/)>2.4 rad (Sec. VII Fl) 


0.173 ±0.015 


eji 




0.170 ±0.015 


eji 


j3 = +0.05 (Sec. VII F 2) 


0.164 ±0.011 


eji 


/5 = -0.05 (Sec. VII F 2) 


0.165 ±0.011 
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VIII. CONCLUSIONS 



Using samples of //// and e/i pairs collected with the CDF experiment during the 1992 — 
1995 run of the Tevatron collider, we have performed a high precision measurement of Xi 
the time integrated mixing probability of 6-flavored hadrons produced at the Tevatron. Our 
measurement, x = 0.152 ± 0.007 (stat.) ± 0.011 (syst.), confirms the trend of all previous 
results from pp colliders, and is significantly larger than the world average % = 0.118± 0.005, 
which is dominated by the LEP measurements at the Z-pole. 
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